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Analyzing the orbital structure and lattice distortions in the CE phase of half-doped manganites, 
we demonstrate that the usual approach directly relating the orbital occupation of Jahn- Teller 
ions to the displacements of neighboring ligands may be misleading. For the correct identification of 
orbital structure, it is necessary to take self-consistently into account the electron-lattice interactions, 
kinetic energy of charge carriers, and crystal-field effects. In certain situations, e.g. in the CE phase 
of single-layered manganite Lao.5Sri.5Mn04, the type of orbital ordering strongly deviates from 
that, which one would deduce from the local lattice distortions. 
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I. INTRODUCTION 

An interplay between spin, charge, orbital, and lattice 
degrees of freedom play an important role in the physics 
of transition-metal oxides, especially those with Jahn- 
Tellcr (JT) ions, such as Mn^+ or Cu^"*". In the materials 
with JT ions an orbital ordering and related lattice dis- 
tortions determine a rich variety of different phenomena, 
and this "orbital physics" attracts now a widespread at- 
tentioni"— . Unfortunately, it is rather difficult to directly 
access orbital state of JT ions. Therefore, the standard 
way to find orbital occupation, widely used for the last 50 
years, is based on structural data, with the assumption 
of the one-to-one correspondence between orbital occu- 
pation of an ion and corresponding JT distortion of the 
surrounding ligands, e.g. Oq octahedral. However, as 
we show below, this straightforward approach sometimes 
fails for the solids with orbital ordering. A well-studied 
layered manganite Lao.5Sri.5Mn04 gives a vivid exam- 
ple of such discrepancy. Below 240 K, this compound 
exhibits a charge ordering characterized by the checker- 
board arrangement of the Mn'^+ and Mn^+ ions. Be- 
low Tjv =110 K, a spin ordering of the CE type appears 
consisting of zigzag ferromagnetic (FM) chains with an 
antiferromagnetic (AFM) stacking between them. The 
CE-type AF was ascribed to an orbital order at the for- 
mal Mn^+ sites^i^. At the same time, the detailed type 
of orbital order involving a spatial redistribution of va- 
lence electrons is still controversially discussed in liter- 
ature. Some experiments are more compatible with the 
3.T^ — r^/Sy^— r^ order— i^ originally suggested in Refs.yS 



while the other suggest the x 



'Iv' 



z^ order——. 



The theoretical treatment based only on electronic 
mechanism of orbital ordering — li^ gives a rather good 
description of half-doped manganites, which reproduces 
the original assignment of occupied "ix^ — r'^ jiy^ — r^ 
orbitals. However, to get a better insight into an ac- 
tual physical situation, one should also include explicitly 



electron-lattice interaction as well as crystal-field effects. 
For undoped manganites, an attempt in this direction 
was recently undertaken in Ref. Il4l . In the present pa- 
per, we analyze the interrelation of orbital order, lat- 
tice distortions, and charge disproportionalization for FM 
zigzag chains in CE phase of half-doped manganites, and 
demonstrate that the standard paradigm of one-to-one 
correspondence of orbital ordering and distortion of near- 
est neighbor ligands fails in these cases. Our results may 
be a "warning sign" that similar discrepancy may also 
occur in other situations, and that the standard way to 
obtain orbital occupation from the local lattice distor- 
tions may be inapplicable in some cases. 



II. A MODEL OF ZIGZAG CHAINS IN 
HALF-DOPED MANGANITES 

A. The model Hamiltonian 

Let us consider an orbital structure in CE phase of 
manganites. In CE phase, the system consists of zigzag 
ferromagnetic (FM) chains, and the neighboring chains 
are ordered antiferromagnetically, as it is shown in Fig. [1] 

This structure can occur e.g. in layered manganites 
such as La2-a;Sra;Mn04 when x ~ 1.5, or perovskite man- 
ganites such as Lai_j;Caj;Mn03 near half doping, when 
X « 0.5. In both these cases, a charge ordering (CO) oc- 
curs with decreasing temperature: we have the checker- 
board CO in the basal plane of, formally, Mn-^^ {t2g^e],) 
and Mn'^^ (^25^6^) ions, with Mn'^^ being strong Jahn- 
Teller ions (one electron on a doubly degenerate e^ or- 
bital), in these systems, there occurs an orbital ordering 
(00) simultaneously with the charge ordering, and the 
00 originally proposed in Refs. [^iGj is shown in Fig. [1] 
As mentioned in the Introduction, it is just this detailed 
type of 00, which is now put under questiorti^, and this 



is the main problem we address theoretieally in our work. 
The model Hamiltonian describing this situation ean 
be written as 
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In this expression, ajj^^ and a„a are creation and an- 
nihilation operators for an Cg electron in orbital state a 
[x^ — y'^ or 2z^ — x^ — j/^) at a site n, Q2n and Q^n de- 
scribe Jahn- Teller (JT) distortions of MnOg octahedron 
(see, Fig. [2|), e^ is the zz component of the tensor, and 
T^'^ are pseudospin operators 
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describing orbital occupation of Cg orbitals; here cr^n''^ 
are Pauli matrices. Note, that in our model the spin of 
eg electrons is always parallel to the core spin S of <2g 
electrons. Moreover, all spins in the chain arc assumed 
to be ferromagnetically ordered, and we omit spin indices 
in Hamiltonian ([1]) . We consider below one zigzag chain: 
as is common in such cases^^ for the CE magnetic order- 
ing of Fig. [T] there will be no electron hopping between 
zigzags with antiparallel spins, and the elastic coupling 
between zigzags, which may be present, does not modify 
the results. 

The first term in Eq. ([T]) is the kinetic energy of eg 
electrons. The hopping amplitudes depend both on the 




Figure 1: (Color online) Orbital and magnetic structures of 
Mn ions in one of Mn02 planes of the CE phase, originally 
proposed in Refs. [^la. A "building block" of the zigzag chain 
consisting of four sites is marked by a dashed line. 




Figure 2: (Color online) (32-type (left panel) and Q3- 
(right panel) distortions of a MnOe octahedron. 
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positions of neighboring sites n and m, and on the orbital 
states of the eg electron before and after hopping. They 
can be written as 
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where minus (plus) sign corresponds to the hopping in x 
{y) direction, see Fig. [1] From the relation between f^ 
for different a and /3, it can be found that the kinetic 



energy term favors the 2x 



and 2y^ 



orbital states in the bridge sites of the CE chain (sites 2 
and 4 in Fig. [T]), while at corner sites (sites 1 and 3 in 
Fig. [T]) the x^ — y^ states turn out to be more favorablei^. 

The second term in Eq. ([T]) is the electron-lattice in- 
teraction, while the third term is the elastic energy. The 
Qzn mode describes the stretching if Qzn > and com- 
pression in the opposite case) of MnOe octahedron in 
z direction, alongside with the compression in x and y 
directions, while Q2n mode describes the stretching if 
Q2n > in a; direction (and compression in y direction 
of MnOg octahedron (see Fig. [2). In both cases, the vol- 
ume of the octahedron does not change. Note that we 
neglect the breathing mode (5i„ in Eq. ([1]), describing 
the uniform stretching or compression. 

In addition to the usual terms describing the electron 
hopping and JT interaction, we introduced in Hamilto- 
nian ([1]) two extra terms, the 4th and the 5th terms in 
Eq. ([T]). £2 describes an external stress acting on MnOg 
octahedra, which tends to elongate them in z direction 
(for Sz > 0) or compress them for e^ < 0. Such term may 
appear just due to crystal structure, or e.g. for films it 
can describe the strain caused by the lattice mismatch 
with the substrate. Specifically, in the layered mangan- 
ites like La2_xSra;M04 there exists the tendency to local 
elongation of MnOe octahedra in z direction (i.e. Ez > 0), 
determined just by the crystal structure and leading to 
such elongation even in the absence of orbital degener- 
acy, e.g. for M=Ni2+ in La2Ni04 (see Ref. [ll). The 
Jahn- Teller effect for ions M with orbital degeneracy will 
operate on the background of this strain due to a layered 
structure. 

The last term in Hamiltonian ([T]) describes an eventual 
contribution to crystal field splitting due to further near- 
est neighbors. Again, for layered materials like A2MO4 



the interaction with further neighbors reduces local sym- 
metry from cubic to (at least) tetragonal, which would 
lead to a crystal-field splitting of eg states even for reg- 
ular MnOg octahedra. In what follows, we will sepa- 
rately study the possible role of strain effects {Ez terms 
in Hamiltonian ([1])) and the longer-range contribution to 
crystal field (A-tcrm in Eq. ([T|)) on the orbital occupa- 
tion vs lattice distortion. 



B. Fully localized case 



When A = 0, E does not depend on 0„, and all states 
\9n) are equivalent. The non-zero A breaks this symme- 
try, making the state \2z'^—x'^—y^) {On ~ 0) when A > 0, 
or the state \x^ — y^) [On = tt) when A < 0, more favor- 
able in energy. In addition, the non-zero stress e^ leads 
to non-zero distortion Qsn even without eg electrons. In 

equilibrium, this distortion is Qa — £z/K ■ The non- 
zero Qa leads to the additional splitting of the eg levels 
with the splitting energy Aq = 'i-gQ^, = 2gez/K. 



First, let us find the orbital structure of the zigzag 
chain in the absence of electron hopping between lat- 
tice sites. A particular orbital occupation for Cg electron 
at site n is caused by the splitting of the doubly de- 
generate eg levels, connected in particular with the JT 
distortions, Q2n and Q^n- The values of Q2n and Qsn 
depend on the orbital state of eg electron at site n. If 
we neglect the intersite hopping, these distortions can 
be found analytically. In this case, the eg electron at a 
site n can be described by the wave function \dn), which 
is the superposition of basic states, \a) ~ \x^ — y^) and 
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Similarly, one can introduce an angle 6' characterizing 
local distortion of MOg octahedra in {Q2, Q3) plane 



tan0' = 



Q2,n 
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In the usual approach, one always assumes that the 
orbital mixing angle is equal to the distortion mixing 
angle 9' , 9 ^ 9' . If so, one could indeed obtain orbital oc- 
cupation from the measure of local M-0 distances, as is 
usually done. We will show below that, generally speak- 
ing, these mixing angles may be different, 9 y^ 9' , i.e. the 
orbital ordering does not necessarily follows local distor- 
tion. 

The energy corresponding to the state |0„) is found 
from Eq. ((!]) (with t"^ = 0) using obvious relations 
Tnl^n) = I - 9n) and T^\9n) = k - 6'„). Minimizing 
the total energy with respect to Q2n and Qsn, we find 
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Already from here, we see that for £2 ^ local distor- 
tions Q2 and Q3 do not exactly follow the orbital occu- 
pation (orbital mixing angle is 6'„, but Q2/Q3 7^ tan0„, 
as it would if the angles 9„ and 9[^ would be the same). 

The total energy then reads 



C. Non-zero hopping: the mean- field 
approximation 

For non-zero hopping, the wave functions |0„) are not 
eigenfunctions of the Hamiltonian ((T|) , and it is necessary 
to diagonalize H numerically or using some approxima- 
tions. In a mean field (MF) approximation, the JT dis- 
tortions are found by minimization of Hamiltonian ([T|) 
averaged over electronic degrees of freedom. It gives 
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This expression is a generalization of Eq. ^ to the case 
of nonzero intersite hopping. As we will show below, the 
orbital occupation characterized by average pseudospins 
(r^) and (r^) is determined not only by the local dis- 
tortion of the first coordination sphere (here, Og octa- 
hedron) but also by other factors such as a contribution 
of further neighbors to the crystal-field splitting, or by 
the effect of band formation. However, given the orbital 
occupation (values of (r^) and (r^)), we can find local 
distortions (the values of Q2 and Q^) from the expres- 
sion (IH), i.e. these local distortions are determined by 
the orbital occupation and by the "external" stress £z- 

Substituting Eq. ([8]) into Hamiltonian ([1]), we obtain 
an effective electron Hamiltonian, which is quadratic in 
electronic operators, and which can be easily diagonal- 
ized. The mean values (r^) and (r^) are found using 
this effective Hamiltonian, and as a result, we find a self- 
consistent equations for (t^'^) and Q2,3n- 

In this procedure, we assume the existence of some 
superstructure in the system. In the case of CE phase, we 
have a one-dimensional effective electronic Hamiltonian. 
The unit cell of the zigzag FM chain consists of 4 Mn 
ions (we enumerate them by Latin subscripts i, j, . . . = 
I, 2, 3, 4, see Fig. [T]), whereas (r^'^) and Q2,3n are equal 
for sites in equivalent positions. However, due to internal 
symmetry of CE chain, the number of sites in the unit cell 
can be reduced to two. Indeed, sites 3 and 4 can be made 
equivalent to sites 1 and 2, respectively, through the use 
of a certain transformation. Note that the JT distortions, 
Q2i and Q2i, obey the following relationships 
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Then, wc introduce new electron operators Cnja 
ing to the fohowing formulas: 
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J = l, 3, (10) 



Thus, we can introduce a new unit cell consisting of two 
Mn ions, and a new electron creation and annihilation 
operators, c„^ and CnA, acting on the states inside the 
two-site unit cell. They can be represented in the form 
of a vector 
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(11) 

where Uj enumerate sites equivalent to the site j, and 
the upper (lower) sign corresponds to j = 1, 2 (j = 3, 4). 
The transformation pO]) for corner sites {j = 1, 3) is 
reduced to the T7i'/2 phase shift in an^a operators, while 
the transformation (fTTj) for bridge sites (j = 2, 4) consists 
of two subsequent transformations: the rotation in the 
pscudospin space by an angle 9 — 27r/3 for j = 2 and 9 = 
— 27r/3 for j = 4, and the additional phase shift. Note, 
that the rotation by an angle 9 = 27r/3 {9 = — 27r/3) 
corresponds to choosing the states |j/^ — z^) and |2a;^ — 
y"^ — z^) {\z'^ ~ x^) and \2y'^ ~ x^ — z^) ) as new basic 
states. After these transformations, sites 3 and 4 become 
equivalent to sites 1 and 2, respectively. 



-^nA 



-lib 



i2b 



) , CnA 




(12) 



where n enumerates now new unit cells. Using relation- 
ships (Uni), (fTTj) . and Eq. ^, we can rewrite Hamilto- 
nian ([T]) in terms of the new electron operators. In the 
momentum representation the effective Hamiltonian can 
be written as 
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and Qsj = Qsi - A/g. 

The spectrum of electrons in CE chain is found by 
diagonalization of this matrix. It consists of four non- 
intersecting bands separated from each other by energy 
gaps. In the case of half-filling, x = 0.5, there is one 
spinless electron for two sites. Since the bands do not 
overlap, only the lowest band is filled at x = 0.5. The 
spectra of electrons at different values of stress Ez are 
shown in Fig. [3l At small Ez, the spectrum of the lowest 
band has a minimum at fc = tt, while at large Ez, k = tt 
corresponds to the maximum. Hence at some value of Ez , 
the spectrum flattens resulting in zero bandwidth. The 
evolution of width Wi of the lowest band with the growth 
of stress is illustrated in Fig. [H 

Let us denote the eigenvalues of e(fc) by Es{k) and its 
eigenvectors by v]l'{k) (s — 1, 2, 3, 4). The mean values 
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where chemical potential fi is found from the condition 
n = l-x = -^(4^c„a)- (16) 



In terms of new electron operators, expressions (l8|) for 
the mean values of JT distortions can be rewritten as 
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Figure 3: (Color online) The spectrum of electrons in CE 
chain es{k) (s = 1, 2, 3, 4 from bottom to top) at small {sz ~ 
0.4, upper panel) and large (ez = 1.5, lower panel) values of 
the stress. Sz is measured in units of Kto/g. The parameters 
of the model are: g'^/Kto = 1.2, K/g = 40, A = 0. fi is the 
chemical potential. In the lower panel, the lowest band s — 1 
is shown in a larger scale. 




Figure 4: (Color online) The width Wi of the lowest band 
(s — 1) as function of strain. The parameters of the model 
are: g^/Kto = 1.2, K/g = 40. 




Figure 5: (Color online) The dependence of 62 (solid red 
curve) and 62 (dashed blue curve) on the stress Sz- The pa- 
rameters of the model are: g^ /Kto — 1.2, K/g = 40. 



III. RESULTS 

Orbital states and local lattice distortions at 
bridge sites 



Eqs. (fT4l) - p8l) together with relationships ([9]) form a 
closed system of equations for finding both JT distor- 
tions and electronic configurations at corner and bridge 
sites of the CE chain. 



First, let us discuss the behavior of orbital states at 
bridge sites. As it will be shown below (subsection D), 
the bridge sites can be characterized by the orbital state 
\9) corresponding to a definite angle in the orbital 
(t'^jT^) plane, whereas for the corner sites, the situation 
is more complicated. Let 62 be the orbital angle for the 
bridge site 2. The corresponding local lattice distortions 
are given by an angle 02 in ((52 1 Q3) plane. The relation- 
ship between 02 and 02 is given by Eq. ([5]), which shows 
that 02 = 02 in the absence of external stress e^ but differ 
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Figure 6: (Color online) The dependence of Q2 (solid red 
curve) and Q'2 (dashed blue curve) on the stress e^ at non- 
zero crystal field splitting A = — O.Sfo- The parameters are 
the same as for Fig. [S] 



Figure 7: (Color online) The dependence of Qi (solid curves) 
and Q'2 (dashed curves) on the hopping integral t at zero and 
non-zero crystal field splitting A. The stress e^ is chosen so 
that gSz/Kto = 0.6. The parameters are the same as for 
Fig. El 



at £2 7^ 0. In Fig. [51 wc demonstrate the plots of O2 and 
6*2 as functions of Sz- At e^ ~ 0, both angles arc equal to 
27r/3 (|2x2 - y2 _ ^2<) orbital). With the growth of e^, 6'^ 
decreases faster than 6*2 . At a certain value of e^ , we have 
^2 = 7r/3. Such 9'2 describes the oxygen octahedron com- 
pressed in y direction and stretched along x and z axes. 
In the usual naive approach with 62 = ^2 , this should give 
\x'^ — z^) orbital lying in the {x, z) plane; this assumption 
was actually used in Ref. [l3 when the authors concluded 
that the orbital occupation in Lao.5Sri.5Mn04 should be 
of x"^ — y^ Iv^ — z^ - type. However, in Fig. |5l we see that 
the the actual orbital angle Q2 at 6*2 = 7r/3 for the chosen 
values of parameters is about 7r/2, that is the electron 
density in z direction is, in fact, much lower than for the 
|x^ — z^) orbital. Physically this is due to the fact that 
the gain in kinetic energy related to the intersite electron 
hopping favors orbitals lying in the (x, y) plane, since the 
electron hopping occurs mainly in this plane. 

Note that by including the term £zQ?,n in Hamilto- 
nian (HI), we took into account a possible contribution 
to a crystal field splitting and, consequently, on the or- 
bital occupation due to distortion of local ligand (here 
Oe ) octahedron. At the same time, the crystal field in- 
cludes also a long-range interaction and it acts directly on 
the splitting of e^ level. This additional splitting, which 
was introduced in Hamiltonian (Jll) in the form — Ar,^, 
favors the orbitals lying in the (x, y) plano^. Indeed, as 
one can easily see, in the layered 214 compounds like 
(LaSr)2Mn04 next nearest neighbors of a Mn ion will be 
four Mn^-^+ in the basal plane (see Fig. |T1), but no such 
ions (or lying at larger distances) in the z-direction. Con- 
sequently, the Coulomb field of these positively-charged 
ions in a6-plane would "pull" the Mn electrons to this 
plane, i.e., they would cause extra crystal field splitting 
stabilizing "flat" x^ — y^ orbitals. 



In Fig. ini we demonstrate the effect of splitting A on 
the behavior of 6*2 and Q'2 as functions of e^. We can see 
that for even relatively small values of A, we have O2 ~ 
27r/3, i.e. the orbital occupation |3y^ — r^ > (or |3a;^ — 
r^ >) even at Q'2 = 7r/3, i.e. for locally compressed (along 
X or y directions) MnOe octahedra. This agrees well 
with the experimental findings and ab initio calculation 
reported in R ef. [g , and explains its discrepancy with the 
results of Ref. [lO- 

As it was already mentioned above, for our layered sys- 
tem the kinetic energy gain determined by the in-plane 
electron hopping integral t favors the orbitals lying in the 
(x, y) plane. In Fig.jTl we show the dependence of 02 and 
6*2 on the hopping integral t. It is clearly seen that the 
difference between 62 and 6*2 grows with t, and 62 — )■ 27r/3 
whereas 62 still remains close to 7r/3 at large t. 



B. Lattice distortions 

Let us now analyze the JT distortions at the corner 
(j = 1, 3) and bridge (j = 2, 4) sites. We will consider 
the case of half-filling, x = 0.5. The numerical calcula- 
tions show that Q21 = 0. In the absence of the breathing 
mode, Qin = 0, we can find the local distortions of Mn 
octahedron in x, y, and z directions. They are related to 
Q2J and Q^j modes by the following formulas 
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Figure 8: (Color online) The local distortions at the corner 
(upper panel) and bridge (lower panel) sites as a function of 
Ez (in units of Kto/g). The parameters are the same as for 
Fig. [5] For corner sites (upper panel) the distortions along x 
and y axes coincide. 



The values of d^'^'^ = 1 correspond to an undistorted 
octahedron. For j = 3, 4, one needs to use relation- 
ships (O- The dependence of the local distortions on the 
stress Ez is shown in Fig. [S] At Sz = 0, for the corner 
site 1 we have compressed octahedron along z axis and 
stretched equally along x and y axes. With the growth 
of the stress, the deformation (elongation) in z direction 
increases, and at some critical value of Sz the deforma- 
tions X, y and z change sign, and the octahedron around 
corner site 1 becomes elongated. 

For the bridge site 2, we have compression in y and z 
direction, and stretching in x direction at small Ez ■ With 
increasing e^. the elongation in z direction increases, 
and deformations in x and y directions decrease. Even- 
tually, for very large e^ , we would have stretched octahe- 
dra along the z axis. But for intermediate values of e^. 
wc may have the situation with the distances (Mn0)2 = 
(MnO)j^ > (MnO)a, (or, for the other edge site, (MnO)2 
= (MnO)j; > (MnO)y), i.e. the locally-compressed (along 
X 01 y directions) MnOg octahedra. As we saw above, in 
Fig. [71 the form of occupied orbitals in this case may still 
be {3y^ — r^) (or (3x^ — r^)). Experimentally this indeed 
seems to be the case in Lao.4Sri.5Mn04^. 




Figure 9: (Color online) The disproportionalization Sn in elec- 
tron density between the bridge and corner sites vs. Ez- The 
parameters are the same as for Fig. [S] 



C. Charge disproportionalization 

Let us now consider the electronic degrees of freedom. 
Note first of all that there exists a charge dispropor- 
tionalization between corner and bridge sites. To mea- 
sure this disproportionalization, we introduce the vari- 
able Sn = (c„2aCn2o) + {Cn2b'^ri2b) - n which dcscribcs 
the deviation in the electron density at the bridge site 
from mean value n. The dependence of Sn on Sz at 
n — 1 — X — 0.5 is shown in Fig. [9l We see, that Sn 
is always positive, that is, the electrons more likely oc- 
cupy bridge sites. We also see that the charge-transfer 
Sn is always less that 1, even for the (unphysically) large 
strain, c/ Ref. UA . 



D. Occupation of orbital states 

Let us now analyze in more detail the orbital configura- 
tions at the corner and bridge sites. First of all, until now 
we considered "pure" orbital states, assuming that the 
electrons occupy a particular orbital. However, strictly 
speaking, for non-zero hopping the on-site orbital state 
can not be described by the orbital wave functions \6) of 
the form of Eq. ([4]): band formation can lead to mixing of 
orbital occupation, and, strictly speaking, it is necessary 
to consider a density matrix p^ (j = 1, 2, 3, 4). In the 
basis of operators a„j„, it can be written as 



Pap 



''Up' 



(21) 



where nj corresponds to the corner site if j = 1, 3, or to 
the bridge site when j = 2, 4. The elements of this ma- 
trix can be calculated using relationships between a„.Q 

and Cnja operators, and Eq. p3)) for (c^^c„b)- The JT 
distortions Q2n and Q^n can be expressed in terms of the 
density matrix. On the other hand the given Q2n and 
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IJQ, 



where plus (minus) sign corresponds to the "Pja (pjb)- We 
will consider 6j = pjb as a measurements of deviation of 
the electronic state at a site rij from the pure state \9j). 
In Fig. [TU] the dependence of 6j and Sj on the stress e^ 
is shown both for the corner (j = 1, upper panel) and 
the bridge (j = 2, lower panel) sites. The angle 9i for 
the corner site changes suddenly from the value 0i = to 
6'i = TT at some critical value of the stress e^ , while the 62 
for the bridge site change continuously from 62 = 27r/3 
to 02 = 0. 

We see that for the bridge site the state 16*2 + tt) is 
almost empty (^2 < 10"'^ <C 1) for any value of e^, while 
for the corner sites the probabilities pia and pu turn 
out to be of the same order of magnitude. That is, for 
the bridge sites, we have nearly pure orbital states 1^2,4) j 
thus justifying our previous treatment, in which we pre- 
dominantly considered orbital occupation and distortion 
of the bridge sites. At the same time, for the corner 
sites both orbitals are populated. This corresponds to 
the conclusions reached in Ref. UM 



E. Discussion of the results 



Figure 10: (Color online) The dependence of 9j (red curves) 
and Sj (blue curves) on the stress Sz calculated for the corner 
{j — 1, upper panel) and the bridge {j = 2, lower panel) sites. 
The parameters are the same as for Fig. [5] 



Qsn can not provide unambiguous determination of or- 
bital states, since the orbital states themselves are mixed. 
In the general case, the density matrix has non-zero 
complex conjugate non-diagonal elements. It can be di- 
agonalized using orthogonal matrix of the form 



Re — 



cos : 



sm ■ 



andia 
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(22) 



(23) 



Transformation (|23[) corresponds to rotation in the or- 
bital space by an angle 6. It is easy to show that the 
value n^^ corresponds to the maximum possible diagonal 

element of the set of matrices of the form ICgfP Rg, that 
is, the orbital state \9j) is mostly occupied. But since 

nA 7^ 0, the orthogonal state \9j + tt) is also occupied. 
The probabilities pja and pjj, to find an electron at a 
site Hj in states \dj) and \6j + tt), respectively, are the 



In our calculations, we modeled the real situation met 
in layered manganites La2-a;Sra;Mn04, especially for the 
composition x = 1.5, for which the system exibits charge, 
orbital, and magnetic ordering of the CE type. We have 
shown that the orbital occupation of the bridge sites 
"Mn'^+" can be treated as an occupation of a particular 
"pure" orbital state, characterized by the mixing angle 9, 
Eq. ([4]), whereas band effects make orbital occupation of 
corner sites less well defined: both orbitals have compara- 
ble occupation at corner sites. This follows just from the 
geometry of the system, with its ferromagnetic zigzags: 
as argued in Ref. Il2l for bridge sites electrons from only 
one orbital can hop to its corner neighbors, whereas at 
corner sites, both orbitals, e.g. \x^ — j/^) and \z^) partic- 
ipate in the band formation. We have also shown that 
the degree of charge transfer in the charge-ordered state 
is (much) less than 1, typically Sn is about 0.2. 

However, the most important is the conclusion that in 
a general case the type of occupied orbitals may strongly 
differ from what one would deduce from the local dis- 
tortion of MOg octahedra. That is, the usually assumed 
one-to-one correspondence between the orbital occupa- 
tion and the local JT distortion, very often used to de- 
termine orbital occupation from structural data, may in 
general break down. 

Two main factors can lead to this effect. One is the 
contribution of further neighbors to a crystal field split- 
ting. Thus, as explained above, in layered materials 



like La2-£i:Sra:Mn04 there are more positively-charged 
Mn'^"'"/Mn "*" ions in xy plane than in z direction, which 
would move the x^ — y^ level down, even though local 
MnOg octahedra may be somewhat elongated in z direc- 
tion, e.g. due to strain, which is always present in layered 
systems. Actually, similar effect was already noticed for 
some other systems, e.g. for BaCoO^i^ or for LaTiO^^. 

Another very important factor is the role of electron 
kinetic energy. For layered systems one would gain max- 
imum kinetic energy if the relevant orbitals are "put in 
plain" - again despite the fact that local octahedra may 
be elongated in perpendicular direction. This factor is 
more important for larger hopping, notably in systems 
approaching localized-itinerant crossover (Mott transi- 
tion). Apparently the considered system La2-KSr2:Mn04 
is in this regime. 

In effect, both these factors lead to the situation, in 
which the orbital occupation may strongly differ from 
that which one would deduce from the structural data. 
Apparently this is what happens in Lao.5Sri.5Mn04, in 
which the MnOg octahedra are distorted so that the Mn- 
O distance to the apex oxygen is large (practically equal 
to the long Mn-0 distance in plane). This distortion 
would in the usual picture lead to the occupation of 
x^ — z^ or y^ — z^ orbitals, which was indeed proposed in 
Ref. [lO- However due to the factors discussed above the 
actual orbital occupation may be quite different, and in- 
deed it was observed experimentally and confirmed by ah 
initio band structure calculations that the occupied or- 
bitals in this case are rather of 2x^— y^ — z^ or 2y^ — x^ — z^ 
type^. Our model calculations show that it is indeed pos- 
sible due to these two factors mentioned above. Thus, 



our results demonstrate that in general the orbital occu- 
pation may strongly differ from that expected from local 
JT distortions - especially in anisotropic systems and sys- 
tems close to the localized-itinerant crossover. Thus, one 
has to be very careful in using the conventional method 
to determine orbital structure from lattice distortions: 
in some cases this often applied method can give quite 
wrong results. 



IV. CONCLUSION 

Through the analysis of the lattice distortions and 
orbital structure of half-doped manganites, we demon- 
strated that in systems with orbital degeneracy and or- 
bital ordering the type of orbital occupation may strongly 
deviate from the one which would be deduced from the 
local JT distortion, i.e. the conventional JT physics may 
be violated in some cases. Thus, one has to be very 
careful in using the standard, widely used method of de- 
termining orbital occupation from the corresponding JT 
distortion. To obtain correct results, it is necessary to 
consider self-consistently orbital occupation, lattice dis- 
tortions, and charge disproportionalization. 
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